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We propose optical phonon lasing for a double quantum dot (DQD) fabricated in a semiconductor 
substrate. We show that the DQD is weakly coupled to only two LO phonon modes that act as a 
natural cavity. The lasing occurs for pumping the DQD via electronic tunneling at rates much higher 
than the phonon decay rate, whereas an antibunching of phonon emission is observed in the opposite 
regime of slow tunneling. Both effects disappear with an effective thermalization induced by the 
Franck-Condon effect in a DQD fabricated in a carbon nanotube with a strong electron-phonon 
coupling. 
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Electrically tunable two-level systems are ideal can- 
didates to study the interaction between fermions and 
bosons under nonequilibrium conditions. Single-qubit 
lasers or nanomechanical resonators [H-Q] are examples 
where concepts from quantum optics, such as the micro- 
laser Q, have been successfully transferred to and com- 
bined with artificial solid-state architectures. Semicon- 
ductor double quantum dots (DQDs) play a similar role 
as model systems with the coupling between electrons 
and the surrounding substrate leading to, e.g., tunable 
spontaneous phonon emission and Dicke-type interfer- 
ence effects [a-llli. 

In this Letter, we propose optical phonon lasing in a 
DQD without the requirement of an additional cavity or 
resonator. We start from the observation that a DQD 
effectively couples to only two LO phonon modes that 
work as a natural cavity. The pumping to the upper level 
is realized by an electric current through the DQD under 
a finite bias. The amplified LO phonons occasionally 
escape from this cavity by decaying into the so-called 
"daughter phonons" [12| that can be observed externally. 
The phonon lasing is possible when the pumping rate is 
much higher than the phonon decay rate r p h. 

We also observe the phonon antibunching in the same 
system when the pumping rate is lower than F p h. We 
emphasize that the phonon statistics can be changed by 
electrically tuning the tunnel coupling between DQDs 
and leads Note that LO-phonon-assisted transport 
through a DQD was theoretically studied by Gnodtke et 
al. |14| and has recently been observed by Amaha and 
Ono [15|. We also note that phonon lasing b y op tical 
pumping was suggested in single quantum dots |l6[. 

Both the phonon lasing and antibunching are weak 
coupling effects that are spoilt by phonon thermaliza- 
tion via the Franck-Condon effect [17Hl9j for a strong 
electron-phonon coupling. In electric transport, the num- 
ber of electrons in the DQD fluctuates, which is accom- 
panied by lattice distortions and thus the creation of 
bunched phonons. We show that this effect is negligible 
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FIG. 1. (Color online), (a) Model for a double quantum dot 
(DQD) coupled to LO phonons. A large bias is applied be- 
tween external leads. The spacing A between the energy levels 
in dots L and R is electrically tunable. When A matches an 
integer (u) multiple of the phonon energy hujp^, the electronic 
state \L) with n phonons is coherently coupled to \R) with 
(n + v) phonons. (b) Phonon mode functions u(r) along a 
line through the centers of quantum dots located at x = ±1Z, 
when the wavefunctions \L) and \R) are spherical with radius 
7Z. The ^-component of u(x, 0, 0) is shown for S(.A)-phonons 
that couple (anti-)symmetrically to the DQD. Note that u x 
is an odd (even) function of x for S (A)- phonons since the 
induced charge is proportional to V • u. 



in DQDs fabricated on GaAs substrates (weak coupling 
case) but surpasses the lasing and antibunching in DQDs 
on carbon nanotubes (CNTs; strong coupling case) [20(. 



Figure QJa) depicts our model of a DQD embedded in 
a substrate, in which two single-level quantum dots, L 
and R, are connected by tunnel coupling Vc- The level 
spacing A between the dots is assumed to be tunable, and 
the total number of electrons in the DQD is restricted 
to one or zero by the Coulomb blockade. The electron 
couples to LO phonons of energy hu> p h in the substrate 
via the Frohlich interaction. Our system Hamiltonian is 
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TL — TL e + "Hph + Hop, 

n e = j(n L -n R ) + V c (dld R + d^d L ), (1) 

n ph = tiw ph J2 N i' ( 2 ) 

?4 P = ^ y^M a|g (og + at q )n a , (3) 

a=L,R q 

using creation (annihilation) operators d^ a (d a ) for an 
electron in dot a and a q (a q ) for a phonon with wave 
vector q. n a = d^da and N q = a q a q are the num- 
ber operators. The spin index is omitted for elec- 
trons. The coupling constant is given by M a>q = 

sj^^il^o) ~ 7@))^( a \ eiq " r \ a )> where I") is the elec " 
tron wavefunction in dot a, e(oo) [e(0)] is the dielectric 
constant at a high [low] frequency, and V is the sub- 
strate volume. The LO phonons only around the T point 
(1 9 1 1/^- with dot radius ft) are coupled to the DQD 
because of an oscillating factor in M a ^ q . This fact justifies 
the dispcrsionless phonons in Hph- We assume equivalent 
quantum dots L and R, whence Mr i0 = Mi J ^ q e lq " rLR 
with r^ R being a vector joining their centers. 

In Hep, an electron in dot a couples to a single mode 
of a phonon described by 



(Eq\M a ,q\ 2 )^ 



(4) 



We perform a unitary transformation for phonons from 
a q to collective phonon coordinates, 



as 



a-L + a-R 
v/20 + 5)' 



a A 



ah - an 
^2(1-5)' 



(5) 



and other modes orthogonal to ag and a a, where S is 
the overlap integral between cll and a R phonons in Eq. 
d4|. Disregarding the modes decoupled from the DQD, 
we obtain the effective Hamiltonian 



H = H e +fiUph N s + \s(as + al)(n L +n R ) 
+fiwph N A + \a{o.a + a\){n L - n R ) 



(6) 



with Ns = agCLSi Na = a A a,A, and dimensionlcss cou- 
pling constants X S /A = (Eq |ML,q±M J?i q| 2 ) 1 / 2 /(2^ph)- 

The mode functions for S- and ^4-phonons are depicted 
in Fig. QJb) along a line through the centers of the quan- 
tum dots. Since the phonons are dispersionless, they do 
not diffuse and act as a cavity including the DQD 21 



A-phonons play a crucial role in phonon-assisted tun- 
neling and phonon lasing, as discussed below, whereas 
5-phonons do not since they couple to the total number 
of electrons in the DQD, til + n R . Both phonons are 
relevant to the Franck- Condon effect. 

Our Hamiltonian H in Eq. © is applicable to DQDs 
fabricated on a semiconductor substrate, where Sw p h = 



36 meV and X s ,a = 0.1 - 0.01 for ft = 10 - 100 nm 
in GaAs [22j. It also describes a DQD in a suspended 



CNT when an electron couples to a vibron, which is a 
longitudinal stretching mode with fta; p h ~ 1 meV, Xa ^ 
1, and As = in experimental situations fill [l9|. 

The DQD is connected to external leads in series, 
which enables electronic pumping. Under a large bias, 
an electron tunnels into dot L from the left lead with 
tunneling rate Tl and tunnels out from dot R to the 
right lead with T R |23[. We also introduce the phonon 
decay rate T p h to take into account the natural decay of 
LO phonons due to lattice anharmonicity [l2[- We de- 
scribe the dynamics of the DQD-phonon density matrix 
p using the Markovian master equation 



■\H, p] + C e p + Cphp, 



(7) 



where C e p = (T L V[d' L ] + T R V[d R ])p and C ph p = 
Tp h (V[a s ]+V[a A ])p, with V\A\p = ApA^ - \{p,A*A} 
being a Lindblad dissipator. Here, we assumed that the 
temperature of the substrate, k R T, is much lower than 

fiWph. 

In the following, we apply the Born-Markov-Secular 
approximation to Eq. ([7]) by diagonalizing the Hamilto- 
nian H and setting up the corresponding rate equation 
in the energy eigenbasis, 
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(8) 



for the probabilities Pi to find the system in an eigenstate 
withLy = (i\[(£ e + £ ph )\j)(]\]\i). The solution of Eq. 
flSJ) with Pi = determines the steady state properties. 

First, we discuss our numerical results in the case of 
Fl,r ~> r p h. We consider A-phonons and disregard S- 
phonons (As = 0). Figure [2ja) shows the current / 
through the DQD as a function of level spacing A, with 
A^ = 0.1 and 1. Beside the main peak at A = 0, we 
observe subpeaks at A ~ vhhjp^ (i> = 1,2,3,.. .) due to 
the phonon-assisted tunneling. At the i^th subpeak, elec- 
tron transport through the DQD is accompanied by the 
emission of v phonons. As a result, the phonon num- 
ber is markedly increased at the subpeaks, as shown in 
Fig.[2{b), in both cases of A^ = 0.1 and 1. However, the 
physics is very different for the two cases, as we will show 
below. 

For Xa = 0.1 and A ~ vfruiph, the electronic state \L) 
with n phonons is coherently coupled to \R) with (n+u) 
phonons 24 , similarly to a microlascr two- level system in 
a photon cavity, if the lattice distortion can be neglected. 
To examine the amplification of A-phonons, we calculate 
the phonon autocorrelation function 



gf{r) = (: N a (0)N a (t) :)/(N A ) 



(9) 



which is the probability of phonon emission at time r 
on the condition that a phonon is emitted at time pHI . 
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FIG. 2. (Color online), (a) Electric current through the DQD, 
(b) ^4-phonon number (Na), and (c) its autocorrelation func- 
tion g A 2 \o) as a function of level spacing A in the DQD. The 
dimensionless electron-phonon couplings are Xa = 0.1 (solid 
lines) or 1 (dotted lines), and As = 0. Io = eTn/(2 + a) is the 
current at A = in the absence of electron-phonon coupling. 
F L = Tb. = 100 r ph and Vc = 0.1 hw ph . 
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FIG. 3. (Color online), (a), (b) Number distribution of A- 
phonons and (c), (d) its autocorrelation function g A 2, '(r) at 
the current main peak and subpeaks (y = 1,3,5). The di- 
mensionless electron-phonon couplings are Xa = 0.1 [1] in 
panels (a) and (c) [(b) and (d)] and As = 0. Note that three 
lines for the current subpeaks are almost overlapped in panel 
(c). r L = T R = 100 T ph and Vc = 0.1 fojph. 



[26j . g A \o) = 1 indicates a Poissonian distribution of 
phonons, which is a criterion of phonon lasing. We thus 
find phonon lasing at the current subpeaks in Fig. [He). 
We mention that these are not changed in the presence 
of finite coupling A5 to S'-phonons. 

When A.4 = 1, the strength of the electron-phonon 

interaction is comparable to the phonon energy. In this 

case, the lattice distortion by the Franck- Condon effect 

severely disturbs the above-mentioned coherent coupling 

between an electron and phonons in the DQD and, as a 

(2) 

result, suppresses the phonon lasing. Indeed, g A (0) > 1 
at the current subpeaks, indicating the phonon bunching. 

To compare the two situations in detail, we show the 
number distribution of vl-phonons in Figs.[J3a) and[3Jb) 
at the current main peak and subpeaks. In the case of 
A.4 =0.1, a Poisson-like distribution emerges at the sub- 
peaks, whereas a Bose distribution with effective tem- 
perature T* is seen at the main peak. T* is determined 
from the number of phonons in the stationary state as 
l^ e hu ph /(k B T*) _ l] = (N A ). When X A = 1, on the other 
hand, the distribution shows an intermediate shape be- 
tween Poissonian and Bose distributions at the subpeaks 
and a Bose distribution at the main peak. In Figs. [3jc) 
and[3fd), we plot the autocorrelation function g A (r) of 
A-phonons as a function of r. In the case of Xa = 0.1, 
g A \r) ~ 1, regardless of the time delay r, supports the 
phonon lasing at the current subpeaks. At the main peak, 
9a\ t ) — 1 + e~ rphr , which is a character of thermal 



phonons with effective temperature T* . When Xa = 1, 
we find an intermediate behavior, g A \r) ~ 1 + S v e~ rphT 
(0 < 8 V < 1), at the ^th subpeak. This indicates that the 
phonons are partly thcrmalizcd by the Franck-Condon ef- 
fect. For larger v values, the distribution is closer to the 
Poissonian with smaller 8 V . 

To elucidate the competition between the phonon las- 
ing and thermalization by the Franck-Condon effect, we 
analyze the rate equation in Eq. |5]). focusing on the cur- 
rent peaks. We introduce polaron states \L(R), n) for an 
electron in dot L (i?) and n phonons with lattice distor- 
tion: 

\L(R),n) = \L(R)) ® T (t) \n) A , T = cr^W), 

(10) 

where T and describe the shift of the equilibrium 
position of the lattice when an electron stays in dots L 
and R, respectively. Note that the lattice distortion pro- 
duces X A extra phonons: (a, ti\Na\o!, n) = n + X A . When 
A ~ vhujph, the eigenstates of Hamiltonian H are given 
by the zero-electron states |0, n) = \0}®\n)A and bonding 
and antibonding states between the polarons 

\±,n) = ^=(\L,n)±\R,n + v)) (11) 

(n = 0,1,2,...) and the polarons localized in dot R, 
\R, n) (n = 0, 1, 2, . . . , v — 1), in a good approximation, 
provided that Vc *C Sw p h- The rate equations for these 
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states are 



Po.n = - r L P , n + ^\A(n\T*\m + v) A \ 2 Pr 



mol.m 



J2 T R \ A (n\T^\m) A \ 2 P R ,, 



m=0 



+ r ph [(n + l)P ,n+i - nP 0t „] 



(12) 



■fmol,n = ^--fmol,n + X] ^ L\A(n\T\m) a\ 2 P(),r 



2 

r p h 



" + 1 + 7T ) Pmal,n+1 ~ ( n +-) P mol, 



(13) 



where P„ 



P+,n + P-,n (n = 0,l,2, 



and 



P R ,n = -r R P fl ,„ + T ph [(n + 1)Pr,„+i - hPh,„] , (14) 

with P RiV = P mo i,o/2 (n = 0, 1, 2, . . . , v— 1). These equa- 
tions yield the current / and the electron number in the 
DQD (n e ) — (til + n R ) in terms of the number of po- 
laroids localized in dot R, (n R ) = Yln=o Pr,u as 



I = er, 



1 



\ n R) 



= l^M. (15) 



2 + a ' 2+a 
with a = r^/r^. The number of A-phonons is given by 



(N A ) = {v + 2\\)-^ + \\(n e ) 



P h 



(16) 



The first two terms in Eq. (|16|) indicate the emission of v 
phonons by the phonon-assistcd tunneling (from dot L to 
dot R) and creation of 2X\ phonons by lattice distortion 
(with two tunnelings between the DQD and leads) per 
transfer of a single electron through the DQD. The last 
term describes the average number of polarons (n e ) in 
the stationary state. 

When T L , R > r ph , wc obtain I = I + 0(T ph /T L . R ), 
where Iq = eY R / (2 + a) is the current at the main peak 
in the absence of electron-phonon interaction, and |27j 



si 2) (o) 



^±^ + o(r ph /r L , R ). (17) 



These explain the numerical results in Fig. [5] at the 
current subpeaks. Equation (fT7f indicates g A (0) ~ 1 
(phonon lasing) for X A <C v and g A \o) — 2 (phonons 
thermalized by lattice distortion) for X\ 3> v. In the lat- 
ter case, the phonons follow the Bose distribution with 
T* for the deduction of (N A ) in Eq. ([T5|> . 

Thus far, we have discussed the phonon lasing in the 
case of Tjj ji 3> T p h. If the tunnel coupling is tuned to 
be Tjj ji < r p h, wc observe another phenomenon, i.e., 
antibunching of LO phonons pjj ]. Figure 2Ja) presents 
a color-scale plot of g A (0) at the first current subpeak 
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FIG. 4. (Color online), (a) Color-scale plot of autocorrelation 
function of A-phonons, g A (0), at the first subpeak of cur- 
rent in a plane of electron-phonon coupling Aa and Tl,r/Tp^. 
T L = Y R = F L , R , As = 0, and V c = 0.1 tkj ph . (b) gf{r) at 
Aa = 0.05 and Yl,r = 0.1 T p h as a function of r (solid line). 

The electric current autocorrelation function S^u'rent ( T ) ls 8 ^ so 
shown by a dotted line. 



in the X A -{T L ^ R /T ph ) plane for r L = T R = T L ^ R and 
X s = 0. At A^ = 0.05 and T LtR /T ph = 0.1, for exam- 

(2) 

pie, g A (0) <C 1, representing a strong antibunching of 
phonons. This is because the phonon emission is regular- 
ized by the electron transport through the DQD. In Fig. 
SJb), we plot the autocorrelation function of the electric 
current 



5currcnt( T ) = ( : n R (0)n R ( T ) :) / '{n R f 



(18) 



where n R is the electron number in dot R. It fulfills 
flcurrent(O) = 0, indicating the antibunching of electron 
transport, since dot R is empty just after the electron 
tunnels out [26|. Remarkably, g A (r) almost coincides 

with .9currcnt( T )- Again at very strong couplings X A > 1, 
neither phonon antibunching nor phonon lasing can be 
observed because of an effective phonon thermalization 
due to the Franck-Condon effect. 

In our calculations, we have neglected the S'-phonon 
coupling that, however, does not affect the dynamics of 
A- phonons as we have checked. We have also disregarded 
acoustic-phonon-assisted tunneling to excited levels in 
dot R because the energy of LA phonons at small wave 
numbers \q\ < 1/7?. is comparable to or lower than the 
spacing between the energy levels for 1Z < 100 nm, and 
the coupling to LA phonons is much weaker than that to 
LO phonons. Indeed, the LO-phonon-assisted transport 
was clearly observed for level spacings A tuned to fia> p h 
and 2ftw p h in recent experiments [l5| . 

Finally, we discuss possible experimental realizations 
to observe LO phonon lasing and antibunching in 
semiconductor-based DQDs. In GaAs, an LO phonon 
around the T point decays into an LO phonon and a TA 
phonon around the L point, which are not coupled to 
the DQD. These daughter phonons can be detected by 
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the transport through another DQD fabricated nearby 
(29l |30| . Alternatively, the modulation of the dielectric 
constant by the phonons could be observed by near-field 
spectroscopy [31j . With a decay rate T ph ~ 0.1 THz in 
GaAs [l2|, however, the lasing condition » r p h 

might be difficult to realize. Other materials with a 
longer lifetime of optical phonons, such as ZnO [32}, may 
be preferable for observing the phonon lasing. 

Our fundamental research of LO phonon statistics is 
also applicable to a freestanding semiconductor mem- 
brane as a phonon cavity [33, 34 1 , in which a resonating 
mode plays a role of LO phonons. Our theory implies 
that a DQD could generate various quantum states of 
mechanical oscillators. 
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In this supplemental material, we derive analytical expressions for the current /, number of phonons {N A ), and 
autocorrelation function of phonons g A (0) at the current subpeaks in Eqs. (15)— (17) in the main material. When 
A ~ vhujph, the energy eigenstates are given by the zero-electron states |0, n) = |0) ® \ri) A , bonding and anti-bonding 
states between polarons |±, n) = ^= (\L, n) ± \R,n + v)) , and polarons localized in dot R, \R, n) (n = 0, 1, 2, . . . , v— 1), 
in a good approximation for Vc <C Sw P h, as mentioned in the main material. We have introduced polaron states 
\L(R),ri) = \L(R)) T w \n) A , with T = e - x A(a\-a A ) _ Tlle density matrix is given by 

GO OO V— 1 

P = Y P o."l°' n ) < ' "I + Y Y P °>"\ a > n > ( a > n \ + Y p R,n\R, n) (R, n\ (19) 

n— cr— ± n— n— 

in the Born-Markov-Secular approximation. We define occupation number operators for zero-electron states, bonding 
and anti-bonding states between the polarons, and polarons localized in dot R as 

OO OO u—1 

no = |0,n)(0,n| = |0)(0|, n mo i = \a, n){a, n\, n R = J2\R,n){R,n\, (20) 

n— <r=i n— n— 

respectively. The relation of no + n m oi + ^i? = 1 holds. The electron number in the DQD is given by n e = n mo \ + tir — 
1 — tiq. The expectation values of the occupation numbers are 

OO OO v — \ 

(n ) = Po,n, (n mo l) = ^ P mo l,n, (™r) = P.R,n- (21) 

n— n— n— 

In the stationary state, the rate equations in Eqs. (12)-(14) in the main material yield 

> p i/ — l 

= -T L Po, n + -y\A(n\T^\m + v) A \ 2 P nlohm + YJ r fl U(n|r t |7n) A | 2 P fl , m + T ph [(n + l)P 0) „+i - nP , n ] , (22) 



2 

m— 771—0 



= -^P mo l,„ + r L U(n|T|77i) A | 2 Po, m + T ph [(«+! + ~) Pnol.n+1 - U + P„ol J , (23) 



2 

m— 

with Pnoi,„ = P+, n + P_, n (n = 0,1,2,.. .), and 

= -T R P Ri „ + T ph [(7i + 1)Pr,„+i - nP R>n ] , (24) 
with P R . V = P mo i, /2 (n = 0, 1,2, . . . , v - l). 

CURRENT AND ELECTRON NUMBER 

First, we express the current / = eFi(no) in the stationary state. For the purpose, we sum up both sides of Eq. 
over n. Using 



VJ \ A (n\T f \m) A \ 2 = A (m\T YJ \n) AA (n\ T^\m) A = 1, (25) 

n=0 V n / 

we obtain 

= -r L (no) + ^(«moi) + T R {h R ). (26) 



With (no) + (n mo i) + (n R ) = 1, we derive 

a 2 

(no) = 77— i 1 + N)= (n-moi) = 77— i 1 ~ 0- + a)(n R ) 
2 + a 2 + a 

with a = Tr/Tl. These equations result in Eq. (15) in the main material: 

1 + (n R ) 2 - a(fi R ) 



I = eT R — -, (7i,, 



2+ a 2+ a 

The summation of Eq. (fiM)) over n yields 

(nR) = -TTj^-Pmolfi- 
21 R 

PHONON NUMBER 

Next, we derive the phonon number which is 



(N A ) = ^ft, n (0,n\N A \0,n) + ^^P a , n (a,n\N A \a,n) + ^PR,v{R,n\NA\R, 



Y nP 0,n + Y ( n+ S + Pmol >« + Y ( n + X a) P ' 



R.n 



I' 

2 

n=0 n=0 n=u 

= (N A n Q ) + (N A n mo {) + (N A n R ). 

We have used 

VN A T = (T f a A T)(T f a A T) = {a\ - \ A ){a A - X A ). 
We multiply both sides of Eqs. (j2"2")) - ([2~i)) by n and sum up over n: 

= -(T L + i»(AUn ) + -y (A^rWd) + Ifl^n^) + ^(tw>, 



= T L (AUn ) - ( -y + r ph J (JYatw) + X A T L (n Q ) + ^ 7^ r « + A l r P hJ (»w) + ^(n*), 

= -(T R + r ph )(N A h R ) + [(y - 1 + Ai)r fl + Air ph ] (n fl ). 
Here, we have used 

f>| A HTV>A| 2 = aHT^ [X |n)AA(n| ) TV)a = Y A(m\TN A -p\m) A . 
From Eqs. (|34[) - (|36p , we obtain Eq. (16) in the main material: 

(N A ) = (v + 2\ A )-^ + \ A (n e ). 
el P h 

PHONON AUTOCORRELATION FUNCTION 

Finally, we derive the phonon autocorrelation function at zero time delay, 

n m (()] _ (--N 2 A :) _ (NX)-(N A ) 
9a W (N A )i (N A )i ' 



where 



v-l 



( n a) = ^Po,n(0,n|iVi|0,n) + ^ P^a, n\N A \a, n) + ^2PR,u{R,n\N%\R,n) (40) 

71=0 er=± n— n— 

oo oo y — \ 

= Yl n2p °-" + £ [ n2 + A ^ 4n + 1 + x a) + v ( n + 5 + 2A a)] p moi,„ + £ h 2 + Ai(4n + 1 + A^)] P R , n 



n=0 



(41) 
(42) 



= (Nino) + (Nln mo i) + (N 2 A fi R ). 
We multiply both sides of Eqs. ((2"2"} - ([2"4"} by n 2 and sum up over n. A similar technique to the last section leads to 



= -(T L + 2r ph )(yV> ) + -f(Nln m oi) + T R (N A n R ) + T ph (N A n ) + ^(N A n mo \) + — 

A -T R + [v + 1 + 8Ai)r ph 



= T L (Nln ) - + 2r ph) (Nln mo i) + 4X 2 A T L (N A n ) + 



v + 4\ 2 



(N A n mo i) 



+ \ 2 A (l + X 2 A )T L (no)+X 2 A 



1 - v - Z\\ 



T R + (1 - v - 6Ai)r ph 



(™mol) - T R (fi R ), 



= -(T R + 2T ph )(Nlh R ) + [4X 2 A T R + (1 + 8A|)r ph ] (N A h R ) 

+ {[(!/- i) 2 + a 2 4 (i - 3\ 2 A )} r R + x 2 A (i- 6X 2 A )r ph } (» fl ). 



(43) 

(44) 
(45) 



From Eqs. (|33 ]) -(|i5 ]l . we find 



(N A )-(N A )=2X A ^(N A n ) + 

1 ph 



t/ + 2A^ r fl , v + <lX 2 A \ , ST N , nx3 [T R 



2 r r 



P h 



P h 

J A (v + QX 2 A 
(n ) - - —{n mol ) 



v + 2X\T L X 2 A (v + QX 2 A ) 



(N A n mol ) + 2X\ + 2j (N A n R ) 

1/(2 - v) Tr 



r p h 



ZX 2 A 



Now we evaluate g A \o) in the case of Tl,r S> T p h. In this case, (tir) = 0(T p h/TL,R.) from Eq. ([2^| . Then 

/ = ^ + o(r ph /r iifl ) 



2 + a 



and 



(N A ) = 



v + 2X\ Tr 



0(1). 



2 + a T ph 
Equations ((33} and ((35} yield 

2(N A n ) - c^JVatw) + 0(1), (N A h R ) = 0(T ph /T L , R ). 
Using (JVa) = (N A n ) + (N A n mo i) + (N A h R ), we have 

(N A n ) = („ + 2X 2 A )-^^ + 0(1), (N A n mol ) = (v + 2A^)- 

(2 + ay r ph 

Using these relations, we obtain Eq. (17) in the main material: 



(2 + a y r ph 



V + ZX A 



3A 2, (0) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



